arXiv:1503.06754vl [nlin.SI] 23 Mar 2015 


On multisoliton solutions of the constant 
astigmatism equation 


Adam Hlavac 

Mathematical Institute in Opava, Silesian University in Opava, Na Rybnfcku 1, 746 01 
Opava, Czech Republic. E-mail: Adam.Hlavac@math.slu.cz 

Abstract. We introduce an algebraic formula producing infinitely many exact solutions 
of the constant astigmatism equation z yy + (l/z) xx + 2 = 0 from a given seed. A construction 
of corresponding surfaces of constant astigmatism is then a matter of routine. As a special 
case, we consider multisoliton solutions of the constant astigmatism equation defined as 
counterparts of famous multisoliton solutions of the sine-Gordon equation. A few particular 
examples are surveyed as well. 


1. Introduction 

In this paper, we continue the investigation of the constant astigmatism equation (CAE) 

z yy + (“) +2 = 0, (1) 

the Gauss equation of surfaces of constant astigmatism immersed in Euclidean space. 
These surfaces are defined as having constant (but nonzero) difference between principal 
radii of curvature. 

Historical roots of the constant astigmatism surfaces can be traced back in works 
of the 19th century mathematicians [4, 5, 6, 7, 21, 17, 16], though the surfaces were 
nameless at that time. For more detailed history we refer to, e.g., [20]. 

In 2009, after about a century in oblivion, the subject of constant astigmatism 
surfaces has been resurrected by Baran and Marvan in the work [3] concerning the 
systematic search for integrable classes of Weingarten surfaces. In the paper, the surfaces 
gained their name and Equation (1) was obtained as well. Recently, the Equation (1) 
has been examined by several authors [13, 19, 14, 18, 15]. 

The main result of this paper is Proposition 1, where an algebraic formula (12) 
provides arbitrary many solutions of Equation (1) from a given seed. The proof is based 
on the observation that superposition formulas [13, Eq. 21] (see also equations (6) below) 
can be conveniently written in a matrix form and, therefore, their iteration coincides 
with the composition of linear transformations, i.e. matrix multiplication. Consequently, 
the same conclusion holds in the case of multisoliton solutions, i.e. solutions of the 
CAE having their counterparts in the well known multisoliton solutions of the sine- 
Gordon equation. Since the algebraic formula for the n-soliton sine-Gordon solution is 
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well known, the corresponding solution of the CAE can be routinely computed using 
Proposition 2. 

The results of Proposition 1 in combination with [13, Prop. 3] also enable 
us to construct arbitrarily many constant astigmatism surfaces by purely algebraic 
manipulations and differentiation once an initial step (including an integration) is 
successfully performed. 

The most important results from the earliest history of the subject of constant 
astigmatism surfaces were, undoubtedly, obtained by Bianchi [4, 5, 6, 7]. He showed 
(see also [21]) that evolutes (focal surfaces) of constant astigmatism surfaces are 
pseudospherical, i.e. with constant negative Gaussian curvature. Conversely, if one 
equips a pseudospherical surface with parabolic geodesic coordinates and takes the 
corresponding involutes, then they are of constant astigmatism. Bianchi also succeeded 
in finding some of the constant astigmatism surfaces explicitly [4, Eq. (30)]. 

A remarkable class of constant astigmatism surfaces was studied by Lipschitz [17] 
and its subclass was later investigated by von Lilienthal [16]. Lipschitz parameterised 
his surfaces by spherical coordinates of the Gaussian image, see Figure 1. Recently, [14], 
we showed that the solutions of the constant astigmatism equation that correspond to 
the Lipschitz class of surfaces, are the Lie symmetry invariant solutions and constitute a 
four-dimensional manifold. The counterpart sine-Gordon solutions are shown to be Lie 
symmetry invariant as well. 



Figure 1 . Lipschitz surfaces of constant astigmatism. 

The aforementioned geometric link to pseudospherical surfaces served as a tool 
for deriving a nonlocal transformation between Equation (1) and the well-known sine- 
Gordon equation 

uj£ V = \ sin 2cu, (2) 

which is the Gauss equation of pseudospherical surfaces parameterised by asymptotic 
coordinates. For detailed description of the transformations see [3]. 

By an orthogonal equiareal pattern (Sadowsky [23, 24]) we mean a parameterisation 
x, y such that the corresponding metric is of the form 

X 

zdx 2 H—dr/ 2 , 
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z being an arbitrary function of x, y. An associated slip line field is a parameterization 
£, rj such that the angle between d x and as well as the angle between d y and d v is 
equal to 7t/4. In [13] we observed that to every surface of constant astigmatism there 
corresponds an orthogonal equiareal pattern on the Gaussian sphere. Note that the 
same result was obtained by Bianchi [7, §375, eq. (20)] in the context of pseudospherical 
congruences. We also showed that solutions of the sine-Gordon equation (2) correspond 
to slip line fields on the same Gaussian sphere. 

The method of [13] for generating solutions of the constant astigmatism equations 
and corresponding surfaces of constant astigmatism has its origin in Backhand 
transformation for the sine-Gordon equation. Let w be a solution of the sine-Gordon 
equation (2). Its Backhand transformation, uA A \ is given by the system 

u;^ A ' ) = uj£ + A sin(u/ A ) + u), co^ = —lo v + — sin(o/ A ) — oo), (3) 

A being called a Bdcklund parameter. The famous Bianchi permutability theorem [5], 
see also [7, 22], and the superposition formula 


tan 


w (AlA2) - u 
2 


Ai + A2 
Ai — A2 


tan 


co 


(Ai) 


CO 


(A2) 


( 4 ) 


enable us to compute solution o/ AlA2 \ the Backhand transformation of u/ Al ) with 
Backlund parameter A 2 , by algebraic manipialations. Therefore, when a seed sohatioai 
co is giveai, the iaitegratioai of the system (3) aieeds to be doaie oarly oaice and the further 
iteration of the Backlund transformation can be perfomred purely algebraically. 

Pairs of siaie-Gordoai solutions related by Backhaaid transformation iaiduce solutions 
of the constant astigmatism equation. I 11 [13, Def. 4] we defmed the associated potentials 
/( A ), a?( A ), o/( A ) correspoaidiaig to a pair of siaie-Gordoai solutions ai,uA A h Iai terms of 
g{\) _ e / (A) |-q ie associated potentials are giveai by 





s (a) a cos (a; W+a;), g^^g^-cos(co^-u), 


\9 W 


sin(u; (A) + co), 


= ^g [x) sin(u; (A) - 00 ), 

A 


A sin(a; (A) + co) _ (A) sin(a; (A) - to) 

? Vr) 


9 


(A) 


A^ A) 


( 5 ) 


Expressiaig = l/g^ 2 iai terms of x^ and y^ oaie obtains a solutioai of the CAE. 

Iai the same paper (Prop. 5) we succeeded iai extending the superposition 
priaiciple (4) to solutions of the CAE. Let co, u/ Al \ co^ X2 \ u/ AlA; d be foiar siaie-Gordoai 
solutions related by the Bianchi superposition priaiciple (4). Theai the associated 
potentials g^ xlX2 \ x^ XlX2 \ y( XlX 0 correspoaidiaig to the pair uh A d,u/ AlA2 ) are related to 
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the associated potentials g^ X2 \ x^ X2 \ y ^ corresponding to the pair u>, by formulas 


gO 1 * 2 ) 
^.(AlA2) 

(AiA 2 ) 


_ ~ A1A2 _„(a 2 ) 

Ai + A 2 — 2A x A 2 cos(u/ Al ) — cu^) 

AiA 2 (a 2 ) 2AiA 2 sin(o; (Ai) — w (A 2) ) (Aa) 

A? - A 2 V Ai + A 2 - 2A x A 2 cos(u/ A b - J ^) 9 

A^~A^ (Aa) _ 2 sin(o/ Al ) — uj ^) 

AiA 2 V g {M) 


( 6 ) 


The above formulas can be regarded as a starting point of our approach in this paper. 

The contents of this paper are as follows. In Section 2 we observe that equations (6) 
can be conveniently exposed in matrix form and the iteration of the superposition 
principle for the CAE reduces to mere matrix multiplication. In Section 3 we handle 
the multisoliton case, i.e. we deal with solutions of the CAE coming from the well 
known n-soliton solutions of the sine-Gordon equation. Section 4 deals with constant 
astigmatism surfaces and slip line fields. In Section 5 we study in detail the n-soliton 
case for n = 1,2, 3 producing exact solutions of the CAE, constant astigmatism surfaces 
and slip line fields, whilst the last section is devoted to the subcase when all Backhand 
parameters are equal to 1. 



( 9 ) 
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The above formulas constitute recurrence relations for the quantities :/f n f ift n \ with 
the initial conditions 

x [1] =x 1: y [1] = y i, g [l] = gi- (10) 


Proposition 1. Let x\,yi,gi be the associated potentials corresponding to the pair 
of sine-Gordon solutions. Let be 4x4 matrices with entries defined by 

formulas 


obi 

D n 


Q b1 

°22 


5. 


bl 

33 


A J+ iAi 

obi _ 
°13 — 

A^ +1 A^ 2sin(u;^ — 

X 2 X 2 ’ 

A j+\ 

Aj_|_i X 1 Aj +1 + A^ — 2Aj + iAi cos(u)^ 

X 2 X 2 

A j+1 A l 

Aj+iAi 

obi 

°24 “ 

—2 sin (a; ^ — ca^), 

1 


—Aj+iAi 


s ti t+i + A? — 2A i+1 Ai cos(«M - 



( 11 ) 


all the other entries being zero. Let 


/ \ 

n ,\n\ 

n— 1 

/ x i\ 

yl J 

"w 

(l/9 W ) 

= Yl s [i] 

i —1 

yi 

9i 

\i /gJ 


( 12 ) 


Then :/f n f if n \ g^ are the associated potentials corresponding to the pair . 

Moreover, if = 1 /, then z^(x^ n \y^) is a solution of the constant astigmatism 
equation (1). 


Proof. One easily observes that relations (9) can be obtained by acting of 4 x 4 matrix 
with the only nonzero entries (11) on the column vector = {x^\y bl, g [j] ,i/g [j] ). 
The recurrent formula (9) with initial conditions (10) then can be conveniently written 
as 

X b+h = 5-blxbI, JW = X 1 = ( Xl ,y u gi , l/ 9l ) T . 

ffence 

n— 1 

X [n] = J-J S [i\ Xi ' 

1=1 


□ 
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x l ‘ = 


71—1 * 

xat ' n 

i=1 Ai +! 




n —1 


-*E 


A i+ iA f- i+1 gW s in(u+ 


n 


Ar>,— 


n—p +1 


i=1 \^i +1 + — 2A^ + iAi cos(a;^ — o;^) p=1 A^_ p+1 A^/ 

l/l TT K+i ~ 


y aF 


n ~ 1 \2 "- 1 ,,[*]'! "-*- 1 \2 _ \2 

fa; — a; J t-t A «-p+i A i 


n 

7=1 


A* 


i+i 


- 2 E 


7=1 


sin 


AP" 1 ^ 


n 

P =1 


-p+1 

An— 


71—p+1 


71—1 


-A, 


■i+1 


[n] _ \n— 1 I I _ 

11 “ = { A 2 +1 + A 2 - 2A i+ iAi cos (a; ^ - a+) 


3. Multisoliton solutions 

Let a+ = 0. Let us define Ajjj := A*, + A;, A^ := Xk — A i and 
a { : = exp( a*£ + ^ + c ») • 

Solving System (3), we get one-soliton solutions 
a;^ = 2 arctan ai, = 2arctana 2 

and, applying the superposition principle (4) to the triple ca+cu+o+l, we easily obtain 
the two-soliton solutions 

[21 o . A| 2 (ai — a 2 ) - [ 2 i 0 , A^ 3 (a 2 — a 3 ) 

a; 1 J = 2 arctan ——-r, o; M = 2 arctan -r. (13) 

Ai 2 (1 + oia 2 ) A 23 (l + o 2 ci3) 

An exact analytic n-soliton solution, in our notation a f n \ of the sine-Gordon equation 
has been obtained by several authors [1, 8, 9, 10, 11, 12], see also [2, 25]. The formula 
best suited for this paper can be found e.g. in [11] and is of the form 

a+] = |arccosp], (14) 

where 

d 2 

(pi = 1 — 2 lndet M (15) 

of ori 

M being the n x n matrix with entries 

Mij — ■ - (\Jo j o :) T / J • 

A i + A j \ Vaittj J 

Note also that oj ^ arises from o+ by increasing all lambdas’ indices by one, namely 
u> ^ = \ arccos (p ^ \ (16) 
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p [n] = 1 - 
M being given by 


d 2 — 

2 In det M 

of oq 


(IT) 


Mij Aj +1 + A 


■3 +1 


Vcti+1%+1 H-7= 

VOjH 


^ i+1®j+1 

Definition 1. By a j-soliton solution of the constant astigmatism equation we shall 
mean a triple (xA'l, gA) formed by associated potentials corresponding to the j- 
soliton solution u)A and the (j — l)-soliton solution (see Diagram (8)) of the 

sinc-Gordon equation. 

Remark 2. To obtain a solution of the CAE explicitly, one would have to express 


^ = l/gli] in terms of xA and yA. However, this is almost never possible in terms of 
elementary functions, see examples in Section 5 and 6. 

A one-soliton solution of the CAE is easy to construct. Following [4, 7], see also 
[13, Prop. 4], x W = x\ and g W = g x can be obtained by differentiation, namely 


9 1 = 


dcu^ 2ai 


dci 1 + a\ ’ 
For y W = y x we have the system 
AiSin^W+o; 101 ) 


[i] _ 


y\ 


V? = 


9 1 

sin(u;W — oh 


Ai gi 

with the general solution 

Tj 

9 i = Aif - — + h, 


xi = 


— Ai, 

1 

“a! 


dln^i a 2 — 1 

dci a, 2 T 1 


(18) 


(19) 


k\ being an arbitrary constant. Setting z\ = 1 fgf, eliminating and dropping the 
lower indices, one reveals the von Lilienthal solution 
1 


z = 


„2 > 


1 — X" 

see Figure 2 on p. 11. 


( 20 ) 


Proposition 2. Let us denote AA = 2 pA ( ph] an y jgA — 2 \J(ipA 2 — 1)(<^^* J — 1), 
where ipA and (pA are defined by (15) and (17) respectively. Then the n-soliton solution 
of the CAE is given by the formula 


f x w \ 

n— 1 

/ -D \ 

yl J 

"w 

(l/9 W ) 

= Yl 

i= 1 

Vi 

9i 

\i /gj 


( 21 ) 
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q\j] 

‘-’ll 


s. 


bl 

22 


5. 


bl 


33 


Aj + i/\i 

^A? 


x 2 

A j +1 


Aj + i - Xj 

Aj+i A] 

1 


bl _ A j+i 


cu\ _ 
°13 — 


A 2 ' 2 






X 2 

A i+i 


Ai + A j+1 

S l f} = -y/2- A® - B® , 

— Aj+iAi 

S$ A 2 + A 2 +1 - A i+1 A! V2 + A® + B® ' 


V 7 2 - A [j] - B [j] 

-X j+1 X xV2 + A^ + B^ 


( 22 ) 


Proof. Formulas (22) follow from plugging (14) and (16) into (11) and employing 
trigonometric identities. □ 


4. Surfaces of constant astigmatism 


Let r be a pseudospherical surface corresponding to a sine-Gordon solution u. Its 
Backlund transformation m A ' is given by the formula 


r (x) = r + 


2A /sin(a; - co (A) ) sin(a; + cu (A) ) \ 
l + A 2 V sin(2w) + sin(2cu) V' 


(23) 


Let us define (cf. (7)) 

r [fc] _ r (AiA2...Afc) j.[fc] _ r (A2A3...Afc + i) 


Then we have the recurrence relation 

2A J+ i (sin(coW -co lJ+1] ) 


.b'+l] - rbl 


= r 121 + 


1 + A i+1 


sin(2o; [jl ) 


„b1 , 


sin(J j] + J j+1] ) 
sin(2o; [2] ) 


.bl 


(24) 


with the initial condition = ro- Surfaces rW are obtained from rW simply by 

increasing all lambdas’ indices by one and replacing uM with u)W. The iteration process 
is shown in the diagram below, cf. (8). 


r 

Ai 

r 


[Q] A LrM Xs > fW A4 >r[3] As > r M — 


Ai 

' 

Ai 

' 

Ai 

Ai 

' 


1] A ' 2 >r 

2] As „ r 

3] A4 „ r 

4] As r 

5] A e 


Substituting A = 1 into (23) one gets what is called a complementary pseudospherical 
surface 


.(i) - 


= r + 


sin w 


sin(2u; 


a/ 1 *) sin^W 1 )) 
r ? H- . — , -r 


sin(2cu) 


v 


Obviously, the surfaces and r^ 2 ^ become complementary when substituting Ai = 1 
into r^. 
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According to [13, Prop. 3], the common involute, r^, of a pair of complementary 
pseudospherical surfaces, r^| Al=1 and is of constant astigmatism and is given by 

the equation 

f[^ = r[l -1 ! — ri^'l lng^|A 1= i, 

where gW is determined by (12) and n^, a unit normal of the constant astigmatism 
surface, is simply 


fjbl _ r bl 


Ai=l 


r b- 1 ]. 


If the surfaces |a 1= i and r m are j-soliton and (j — l)-soliton pseudospherical surfaces 
respectively, then the corresponding common involute, r^, will be called a j-soliton 
surface of constant astigmatism. 

Let us also remark that ri ^ (£, rj) parameterises a unit sphere by slip lines (see the 
Introduction, for details see [13]). 


5. Examples of multisoliton solutions 

In this section we provide explicit formulas for some multisoliton solutions of the 
constant astigmatism equation as well as corresponding constant astigmatism surfaces. 
Firstly, let us introduce a notation. Let us define 

a:=£r~r}, /3:=f + V 

which is nothing but the space and time coordinates in which sine-Gordon equation is 
of the form upp — u Qa = sin a;. Also recall that 

a % := exp^Aj^ + -j- + c^, X kl := A*, + A/, X kl := A& — A; 

and, in order to have short formulas, let us define 

a := exp(£ + 77 + c), A® : Xf. — Xf, X®:=X 2 k -Xf, 

K := ^ +1) K := ^ := := ^l - 1- 


5.1. One-soliton solutions 


In Section 3 we constructed the one soliton solution of the CAE corresponding to the 
pair u = 0 and a;! 1 ] = 2 arc tan a\. It belongs to the von Lilienthal class. Let us proceed 
to the corresponding surfaces of constant astigmatism. The family of well known one- 
soliton Dini’s surfaces is 




AX 2 a2Cosa 4A2a2sina 2X 2 (al — 1)\ 

^ t( a 2 + 1 ) ’ ^2 ( a 2 + 1 )’ -^2 ( a 2 + 1 ) ' 
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r m = 


/4AiOiCosa 4Aiaisino; 2Ai(a 2 — 1) 


VA|(a? + l)’ Afta^ + l) 

we obtain the pseudosphere 




Ai"(oi + 1) 


r [11 U 1=1 = 


2 a cos a 2 a sin a 1 — a" 
2 , i > 2 , i )P 4" A . 2 


a + 1 7 a z + 1 


1 + at 


(25) 


Note that in this case the ‘seed surface’ = r 0 is degenerated and coincides with the 
z-axis (0, 0, (3). 

The Gaussian map of corresponding constant astigmatism surface is 


n [i] = r [i )| Al=1 - r 0 = 


2 a cos a 2 a sin a 1 — of 


2 i i ’ 2 i i i i 2 / ’ 

& H - 1 & H - 1 1 ~h d / 

forming the net of 45° loxodromes on the unit sphere. The family of one-soliton constant 
astigmatism surfaces is then 

2 ka 


/ -2a 


■W = r 0 - n [i l ln(/c 5 f [1 ])| Al=1 = 


, In i-o ) cos a 

1 + a z V1 + a 

—2a ( 2 ka 

m 


1 T a 

P + 


l + a z 


sm a 


a 2 - 1 




2/ca 


(26) 


k being a real constant. The surfaces coincide with the von Lilicnthal class [16], see 
Figure 2; for detailed description and pictures see [3]. Evolutes of the surface r 41 are 
the pseudosphere rW| Al=1 and the z-axis r 0 . 


5.2. Two-soliton solutions 

The two-soliton solutions of the sine-Gordon equation are given by (13) and the two 
soliton solution of the CAE corresponding to the pair and is 


x [ 2 ] _ A2A1 _ A^ (o| ~ o|) + A 12 (a 2 a2 — 1 ) _ 

A12 A^2 (a 2 + a^) + A 12 (a 2 a 2 + 1 ) — 8AiA2aia2 


2 ^„2„2 


y p] = 


A® 


12 


(Ak - h) + 


2(1 + aia 2 )(ai — 02 ) 


AiA 2 v " ' a^l + a 2 ) 
[ 2 ] —2AiA 2 ai(l + a 2 ) 


J = 


z™ = 


A 12 ( a i T a|) T A 12 (a|a <2 + 1) — 8 AiA 2 aia 2 
1 / A 12 (a? + O 2 ) 4” A 12 (a 2 a 2 + l) — 8 AiA 2 aia 2 \ 2 


9 


[2p 


2AiA 2 ai(l T a 


7 


(27) 


The triple zl 2 ], x^ 2 \ y ^ provides a solution of the CAE in parametric form; the plot can 
be seen in Figure 3. Eliminating rj one obtains an implicit formula for the function 
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Figure 2. Pseudosphere rM(£,ry)|^ 1= i (upper left); von Lilienthal surface (£,/?) for 
lnfc = 0.3 (upper right); Gaussian map riW^,??) (lower left); solution z = 1/(1 — x 2 ) 
of the CAE (lower right). 


z(x,y ) = z^ 2 \x^ 2 \y^), namely 
y — 2 In a 2 

where 


® 2 lnai ^ 2(1 + aia 2 )(ai - a 2 ) 


Ai A 2 


Ol(l + ct 2 ) 


a i = 


a-2 — 


-(a; 2 Af 2 — A 2 A 2 ) 2 z 2 — 2A® 2 (x 2 A 12 — A 2 A 2 )^ + 2K AiA 2 A® 2 \fz — Af 2 
(a; A® + AiA 2 ) 2 (4A?A 2 ^i + Kz) + 4A 2 A 2 A® 2 ^ + A'A® 2 


K = 16A 2 A^ - [A® 2 + (a; 2 A® 2 - A 2 A 2 )^] 2 . 


4A 2 AfV- + K 


A® 2 + (z 2 A® 2 - A?A 2 )/ 

Using (24) one obtains the two-soliton pseudospherical surface 


r [ 2 ] = ^12 


\® \® E> 
A^ A 2 ±ii 


'Ri 


A, 


i ?3 sin a + —i? 2 cos a 


0 


0 


/ 


+ 


0 

0 


\ 


Ra 


\ 0 + \f\fRj 
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Figure 3. Two soliton solution (top) and three soliton solution (bottom) of the CAE, 
Ai = 1.2, A 2 = 1.5, A 3 = 1.8, Q = 0. Right parts of the figures show the behavior around 
the origin. 


where 

R\ = A + 2 (af + o 2 ) + A 12 (01O2 "h l) — 8 AiA 2 a 2 ai, 

i ?2 = 2AiA 2 (l + aia 2 )(ai — af), 

( 28 ) 

R 3 = A 2 A ^ o 2 (l + a 2 ) — AiA®ai(l + a 2 ), 

R 4 = 2Af 2 A+ 2 (AiA2 — l)(a 2 — a 2 ) — 2A®2A 12 (AiA 2 + l)(a 2 a2 — 1)5 
for picture see Figure 4. The Gaussian map of corresponding constant astigmatism 
surface is 


hPl = rM 


f Ni sin a + N 2 cos ob 
a 1= i — r W = | N 2 sin a — N\ cos a | , 

N 3 


where 

„ _ „A?ff' 2 _ „A? R 3 4o 2 A 2 _ K 2A 2 (a 2 - 1) 

1 Afif',’ 2 A®if', A® (a 2 + 1)’ 3 2A®ff; + A®(a 2 + 1) ' 

Obviously, R’ t arises from Ri by substituting Ai = 1, namely 

R i — A^~^(o 2 T o 2 ) T A9 (gCc?A T l) — 8A2UO2) 

R 2 = 2A 2 (1 + aa 2 )(a - a 2 ), 

R's = ~ Af a (l + a l)i 

R'a — — 2Xf\f(X 2 — l)(a 2 — a 2 ) — 2AfA 2 (A 2 + l)(a 2 a 2 — 1). 


( 29 ) 
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Figure 4. From the left: One-, two- and three-soliton pseudospherical surfaces, r W, rt 2 ], rW 
respectively, Ai = 1, A 2 = 1.5, A 3 = 1.8, Cj = 0. 


Recall that rb 2 ](£, rf) parameterises the unit sphere by slip lines, for a picture see the left 
side of Figure 5. One can easily check that in the case when q = 0, the slip lines net is 
symmetric with respect to the point h^(0,0) = (1,0,0), i.e. if = ( n x ,n y ,n z ), 

then rb 2 l(— £, — 77 ) = (n x ,—n y ,—n z ). In the picture, we observe the sphere from the 
positive part of the rc-axis and the symmetry of the pattern can be clearly recognized. 
Some parts of the sphere are multiply covered and one can see singularities (envelopes) 
of slip lines. We leave open the question (suggested by a few numerical calculations) 
whether coordinate lines f = const and g = const converge to the poles, i.e. points 
( 0 , 0 , ± 1 ). 




Figure 5. Slip line field fit 2 ], A 2 = 1.5, q = 0 (left) and its limit n^ 2 }, c = 0 (right) with 
coordinate lines £ = 0 and g = 0 highlighted (thick black curves). 


To obtain an associated orthogonal equiareal pattern one needs to invert the 
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transformation (x,y) -B- (£, 77 ) given by first two equations in (27), which is not possible 
in terms of elementary functions. The parameterisation rd 2 l(a;, y) then would provide 
orthogonal equiareal net sought. 

The corresponding family of two-soliton constant astigmatism surfaces having 
evolutes r| a 1= i and a Dini surface fM is then 


f [ 2 ] = jd 1 ] _ f nPl = 


(R sin a + R 2 coscA 
R 2 sin a — Ri cos a 

V P — R3 j 


where 


~ \©zy ~ 

/ = ln(%l 2 l)| Jl=1 , fli =2 = 2 / 

5 _ 1 , #4 . 2V(a 2 -l) 

R3 ~2 i ^W l + Xf(4 + 1) u + 1) - 


X, R, 


+ 


4 A 2 C 12 


A®i?) Af(a 2 + 1) 


(/ + I); 


A part of the surface can be seen in the left side of Figure 6 . One can observe cuspidal 
edges obviously related to singularities of corresponding slip line field (Figure 5). 



Figure 6. Two soliton surface ft 2 ] (left) and three soliton surface ft 3 ] (right) of constant 
astigmatism, A 2 = 1.5, A 3 = 1.8, c* = 0, k = 1. 


5.3. Three-soliton solutions 

The 3-soliton solution, y^\ f/ 3 !), of the CAE corresponding to the pair 

[3] _ 2 arctan ( ^12-^13-^2 3ai ~ ^lAl3^23 a 2 + '^ 12 '^ 13*^23^3 + A 12 A 13 AggOjOgOa 
' / ^12'^13'^23 a l a 2 — ^I2^13-^23 a l a 3 T -^I2'^13'^23 a 2 a 3 + ^12^13^23 
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and can be easily obtained using Proposition 2, explicit formulas being too lengthy 
to be written here. The graph of the solution z (x,y) = z^(x^,y^) can be seen in 
Figure 3, a multivaluedness of the function z being clearly identified. In the right side of 
the figure (values of z near the point (0, 0, 0)) one can observe that at least eight values 
of z may correspond to one particular choice of x and y. 

Proceeding to the constant astigmatism surface we firstly compute the three-soliton 
pseudospherical surface (see Figure 4) 


,[3] = r [2] + 


2A 3 /sin(a;^ — u;^) r 2 i sin(o;^ + 


A? + 1 


sin 2a; ^ 


r[ 2] + 


sin 2a; ^ 


,[ 2 ] 


The Gaussian map of corresponding constant astigmatism surface is then 


fj[3] = r [3] 


Ai=l 


_=[ 2 ] 


forming a net of slip lines, for picture see the left side of Figure 7. 



Figure 7. Slip line field nl 3 ], A 2 = 1.5, A 3 = 1.8, c, = 0 (left) and its limit rii 3 }, c = 0 
(right) with coordinate lines £ = 0 and y = 0 highlighted (thick black curves). 


Finally, we construct the family of constant astigmatism surfaces having evolutes 
G 3 1 | a 1= i and r^. They are 

ft 3 ! = rf 2 ] — fiI 3 ! ln(fc^)| Al=1 , 

k being a constant. A picture of rather complicated surface can be seen in the right side 
of Figure 6. 

6 . Examples of multisoliton solutions with A = 1 

Computing the limit for Aj —> 1 allows us to iterate the Backlund transformation 
with A = 1 and construct solutions of the sine-Gordon equation and corresponding 
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pseudospherical surfaces. Slip-line fields are also available. However, taking the limit for 
(Ai,A 2 ) —> (1,1) and (ci,c 2 ) —> (c, c) in (27) yields 

lim = oo, lim = 0 , lim < 7 ^ = oo. 


Fortunately, finite limits result when applying the scaling symmetry 


S 8 {x,y,g) ->■ ( sx,y/s,sg) 

of the CAE to the triple (x^, y^ 2 \ g^), the constant s being suitably chosen. In what 
follows we iterate the Backlund transformation with A = 1 (the case examined e.g. in 
[25]) and we find the corresponding solutions of the CAE as the limits of the solutions 
obtained in previous section. 

Let us denote the k -th iteration of the Backlund transformation with A = 1 by uj^ k \ 
according to the diagram 


u 


{ 0 } 


■u 


{ 1 } 


■u 


{ 2 } 


■u 


{3} 


■U 


{4} 


The corresponding pseudospherical surface is, obviously, denoted by r^. We also extend 
the notation to the solutions of the constant astigmatism equation, i.e. x^ k \y^ k \g^ k \ 


6.1. One-soliton solutions 

Starting with the zero solution 0 /°^ = 0, the one-soliton solution of the sine-Gordon 
equation is 

= 2 arctana = 2 arctan e^ +T?+c 

and the one-soliton pseudospherical surface, is the pseudosphere (25). 

The one-soliton solution of the CAE can be easily found by substituting Ai = 1 
into Equations (18) and (19). Actually, since the solution does not depend on Ai (it 
disappears when eliminating parameters rj from (18) and (19)), we obtain precisely 
the von Lilienthal solution (20). 

Corresponding constant astigmatism surfaces are exactly those given by (26). 


6.2. Two-soliton solutions 

Taking the limit of for \ —» 1 and q —>■ c yields 


u { 2 } _ 2 arc tan 


2 aa 


= 2 arctan 


2e «+^( ? - .,) 


1 + a‘ 1 + e 2K+,+c) 

The corresponding two-soliton pseudospherical surface is 


r {2} = 


/ 4aa(l + a 2 ) sin a + 4a(l + a 2 ) cos a 
- 7—7 -- 4a(l + a 2 ) sin a — 4aa(l + a 2 ) cos a 

+ 2(2«W + 1 \p(f + 2 ^ +l y + \)-2(a*-l), 
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fi{2} _ r {2} _ r {i} _ _ 1 _ 

a 6 + (4a 2 + 3) (a 4 + a 2 ) + 1 

/4a(a 2 + l) 2 asina + 2a[a 4 + 2(1 — 2a 2 )a 2 + l] cosa\ 
x I 2a[a 4 + 2(1 — 2a 2 )a 2 + 1] sin a — 4a(a 2 + l) 2 acosa J , 

\ (1 - a 2 ) [a 4 + 2(1 - 2a 2 )a 2 + 1] ) 

see the right side of the Figure 5 . The picture indicates some symmetries of the pattern. 
Indeed, if c = 0 and nl 2 l(£, 77) = ( n x ,n y ,n z ), then 

• n^ 2 ^(—£, — 77) = (n x ,—n y ,—n z ) (symmetry with respect to the point n^^OjO) = 
(1,0,0)) , 

• n {2} (h,0 = n z ) (symmetry with respect to the zero meridian - the image 

of the line £ = 77 ), 

• rC 2 l(—77, — £) = (n x ,n y , —n z ) (symmetry with respect to the equator - the image 
of the line f = —77 ). 

Coordinate lines of the slip line field n^ 2 l converge to the poles, namely 


lim = lim = (0,0, =1=1). 

£—»=too r/— >±00 

According to the definition of slip lines, the corresponding orthogonal equiareal pattern 
has the same symmetries and the same limit behavior. Unfortunately, no picture is 
provided; to do that one has to reparameterise the sphere which requires expressing £, 77 
in terms of x^ 2 \y^ 2 \ i.e. inverting the transformation described by first two equations 
in ( 30 ). Plugging the functions £(x,y),r](x,y) into nl 2 l(£, 77) then would provide the 
orthogonal equiareal net h^(x,y). 

The two soliton solution (x^ 2 \y^ 2 \ g^) of the CAE can be obtained by taking 
a limit for X t —* 1 and q —>■ c of the triple (x^ 2 \ y^ 2 \ g^) scaled by the factor 
s = — (A 2 - A1) 2 , i.e. 


(x { 2 } ,y { 2 } ,g {2} ) = lim 5 _ (A2 _ Al)2 (a; [2l ,j/ [2l ,^ [21 ). 
a i y 1 

Ci—>C 


We have 


4e 2(^+c)(^ - 77) 

e 4 (H+rj+c) + 2 [ 2 (f - 77) 2 + l]e 2( ^ +,?+c - ) + 1’ 

e 2(M-7?+c) ^ + 77 - (£ - 77) 2 + c] + £ + 77 + (£ - 77) 2 + c 

e 2 (£+» 7 -Tc) ^ ’ 

2e Wr?+c(e 2( «+r)+c ) + 

e ^+v+c) + 2 [ 2 {f - -77) 2 + l]e 2< ^ +r?+c ) + 1 ■ 


( 30 ) 
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Then z ^ = 1/g^ 2 as a function of x^ 2 \y^ is a solution of the CAE. Eliminating 
parameters rj provides the implicit formula for the function z = z(x,y) = 
z { 2 } ^ x { 2 }, y{ 2 }^, namely 

If 2 — If xz \ 2 zx 2 + 1 

V= L 2 + 1\ZX 2 + l) _n ’ = y/ Z -{ ZX * + lY -yf Z ' 

Finally, we write down the constant astigmatism surface having evolutes A 1 1 and 
A A For brevity, let us denote ln(fcr/ 2 l) by /. Then 


f {2} = r {i} 


/n ,2 > = 


(R sin a + R-2 cos a\ 
R-2 sin a — Ri cos a 

V P + R 3 


where 

~ 4a(a 2 + 1 )fa 

1 “ _ a 4 + 2(2a 2 + l)a 2 + l’ 

5 _ 0 (/-l)a 4 -2[2a 2 + (2a 2 -l)/ + l]a 2 + /-l 

2 ~ a 6 + (4a 2 + 3) (a 4 + a 2 ) + 1 

~ _ (/ — l)a 6 ~ [(4a: 2 — 1)/ + 4 a 2 + l](a 1 — a 2 ) — / + 1 
3 a 6 + (4a 2 + 3) (a 4 + a 2 ) + 1 


The results from this example coincide with ones obtainable using a reciprocal 
transformation for the constant astigmatism equation introduced in [15], albeit the 
parameterisation of the results is different. For instance, one can observe apparent 
similarity between the surface in Figure 8 and the surface [15, Sect. 8, Ex. 5, 
Fig. 3], 



Figure 8. Two soliton surface, A 2 }, of constant astigmatism, k = 1, c = 0. The central 
part of the surface is zoomed in the right. 
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cn {3} - 2 arctan a ( 2a + W + 1 ) + a' 
u -2arctan a2(2a;2 _ 2 ^ + i) + i 


or, alternatively [25], 

r 3 i _ (1/3 + 1 sinh 2(j3 + c) + cosh 2 (/3 + c) + 

1 = —2 arctan - • —-f —-.-rvw ,— 

\a (3 + \ sinh2(/3 + c) — coslr(/3 + c) — 

The corresponding three-soliton pseudospherical surface is 


a 


a 


jd^} — 


{2} ( smjuW -u^) A2} ^ sin^W 31 ) 


sm 2 cu {2} ^ + sin 2 u; {2} 

and the slip-line field on the constant astigmatism surface’s Gaussian sphere is given by 

fj{3} = r {3} _ r {2} 


T 2 } 


formulas being too lengthy to be written here, hence omitted. However, the picture can 
be seen in the right side of Figure 7 and the properties of the pattern (symmetries when 
c = 0 and limits of the coordinate lines for £, r) —y ±oo) are exactly same as in the two 
soliton case from the previous subsection 6 . 2 . 

In order to obtain the three soliton solution (x^ 3 \y^ 3 \ g^) we rescale the triple 
(xl 3 !,y[' 3 l, r/1 3 ]) by the scaling factor s = (A 2 — Ai) 2 (A 3 — Ai ) 2 and then compute the limit 
for A* —» 1 and c t —> c, i.e. 


y 1 


Ci—>C 


The result is 


is 


2.(3} _ 


y { 3} = 


—4a 2 [(a 2 — f3)a 2 — a 2 — j3\ 

a 6 + [4(a 2 — f3 ) 2 + 8 a 2 + 3]a 4 + [4(a : 2 + f3 ) 2 + 8 a 2 + 3]a 2 + 1 ’ 

2 a (a 2 — 3 f3 + §)a 4 — 2(a 4 — a 2 + 3 (3 2 — |)a 2 + a 2 + 3/3 + | 
1 a 4 + 2(2a 2 + l)a 2 + 1 : 


S {3} = 


2 a [a 4 + 2 ( 2 a 2 + l)a 2 + 1 ] 


a 6 + [4(a 2 — (3) 2 + 8 a 2 + 3]a 4 + [4(a 2 + (3) 2 + 8 a 2 + 3]a 2 + 1 
Finally, the constant astigmatism surface r { 3 }, part of which can be seen in Figure 9, 
f{3} =r {2}_ ln (^{3} )il {3} ) 


k being a real constant. 
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Figure 9. Three soliton surface, r i 3 }, of constant astigmatism, k = 1, c = 0. 
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